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Abstract. As an alternative explanation of the dimming of distant supernovae it 
has recently been advocated that we live in a special place in the Universe near the 
centre of a large spherical void described by a Lemaitre-Tolman-Bondi (LTB) metric. 
In this scenario, the Universe is no longer homogeneous and isotropic, and the apparent 
late time acceleration is actually a consequence of spatial gradients. We propose in 
this paper a new observable, the normalized cosmic shear, written in terms of directly 
observable quantities, and calculable in arbitrary inhomogeneous cosmologies. This will 
allow future surveys to determine whether we live in a homogeneous universe or not. 
In this paper we also update our previous observational constraints from geometrical 
measures of the background cosmology. We include the Union Supernovae data set 
of 307 Type la supernovae, the CMB acoustic scale and the first measurement of the 
radial baryon acoustic oscillation scale. Even though the new data sets arc significantly 
more constraining, LTB models - albeit with slightly larger voids - are still in excellent 
agreement with observations, at x^/d.o.f. ~ 307.7/(310— 4) — 1.005. Together with 
the paper we also publish the updated easyLTB codeQ used for calculating the models 
and for comparing them to the observations. 
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1. Introduction 

The nature of the matter and energy contents of the Universe is still a mystery. While 
many cosmological parameters are known to a few percents precision - which suggests 
we can start to speak of a Standard Cosmological Model - some of these parameters 
correspond to quantities for which we completely ignore its nature. According to this 
model, only 4% of the Universe is made of something whose nature we know. What is 
the rest made off? Moreover, our current model of the Universe is based on a set of 
symmetries we have never tested very thoroughly, like homogeneity and isotropy, and 
some assumptions we have taken for granted, perhaps too lightly, like the generalized 
Copernican Principle, also known as the Cosmological Principle. These symmetries and 
assumptions provide the theoretical basis for the nowadays predominant homogeneous 
and isotropic Friedmann-Robertson Walker (FRW) model of the Universe. 

However, distant supernovae appear dimmer than expected in a matter-dominated 
FRW Universe. The currently favored explanation of this dimming is the late time 
acceleration of the Universe due to a mysterious energy component that acts like a 
repulsive force. Observations seem to suggest that this so-called Dark Energy is similar 
to Einstein's cosmological constant, but there is inconclusive evidence. There has been a 
tremendous effort in the last few years to try to pin down deviations from a cosmological 
constant using deep galaxy catalogues like 2dFGRS [1] and SDSS [2], and extensive 
supernovae surveys like ESSENCE [3], SNLS [1], UNION [3] and SDSS-SN [H], and 
many more are planned for the near future e.g. DES PAU P, [9], BOSS [lO] and 

LSST [n]. 

In the meantime, the realization that the universe around us is actually far 
from homogeneous and isotropic has triggered the study of alternatives to this 
mysterious dark energy. Since the end of the nineties it has been suggested by various 
groups [T2l[l3l[Tl[T5l[T6l[I7l[T8l[T9l[2ni[2ll^ that an isotropic but 

inhomogeneous Lemaitre-Tolman-Bondi (LTB) universe could also induce an apparent 
dimming of the light of distant supernovae, in this case due to local spatial gradients in 
the expansion rate and matter density, rather than due to late time acceleration. It is 
just a matter of interpretation which mechanism is responsible for the dimming of the 
light we receive from those supernovae. Certainly the homogeneous and isotropic FRW 
model is more appealing from a philosophical point of view (it satisfies the Cosmological 
Principle, i.e. it assumes spatial sections are maximally symmetric), but so was the static 
universe and we had to abandon it when the recession of galaxies was discovered at the 
beginning of last century. 

There is nothing wrong or inconsistent with the possibility that we live close to 
the centre of a gigaparsec-scale void. Such a void may indeed have been observed as 
the CMB cold spot [271 [28t [29] and smaller voids have been seen in the local galaxy 
distribution EI]. The size and depth of the distant voids, i.e. tq ~ 2 Gpc and 
VLm ~ 0.2, within a flat Einstein-de Sitter universe, seems to be consistent with that in 
which we may happen to live [32] , and such a local void could account for the supernovae 
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dimming, together with the observed Baryon Acoustic Oscillations (BAO), the CMB 
acoustic scale, the age of the universe, the local rate of expansion, etc. |23j . 

Observations suggest that if there is such a large void, we should live close to the 
centre, otherwise our anisotropic position in the void would be seen as a large dipole in 
the CMB [33]. Of course, we do observe a dipole, but it is normally assumed to be due to 
the combined gravitational pull of the large scale structure in the local universe, such as 
the Virgo cluster, and the Shapley super cluster. There is always the possibility that we 
live off-centre and we are moving towards the centre of the void, so that the two effects 
are partially cancelled, giving rise to the observed dipole. However, such a coincidence 
could not happen for all galaxies in the void and, in general, clusters that are off- 
centred should see, in their frame of reference, a large CMB dipole. Such a dipole would 
manifest itself observationally for us as an apparent kinematic Sunyaev-Zeldovich (kSZ) 
effect for the given cluster. The still preliminary observations of the kSZ effect allowed 
us recently [21] to put some bounds on LTB models. It is expected that in the near 
future, improved measurements of this effect with the Atacama Cosmology Telescope 
(ACT) [31] and the South Pole Telescope (SPT) [3S] will very strongly constrain the 
LTB model. 

In fact, what we require are new data sets that could be used to constrain further 
the LTB models and eventually rule them out, if possible. We thus use the recently 
released compilation of 307 supernovae by the Supernova Cosmology Project - UNION 
collaboration [S]. Moreover, in this paper we make use of a new and very interesting 
observation of the baryon acoustic oscillation along the line of sight, otherwise known as 
radial BAO (RBAO) [HUl [37] . It provides a direct measurement of the rate of expansion 
along the line of sight, Hl{z), whose integral determines the luminosity distance and 
could in principle be different from that transverse to the line of sight, Hxiz), which 
relates to the angular diameter distance. Both appear in the Einstein equations and it is 
possible to disentangle their respective effects by looking at a new observable, the ratio 
of shear to expansion in LTB models. In FRW universes, the shear of the background 
geometry is identically zero, while in LTB models it can be significantly different from 
zero, and becomes maximal at a certain redshift, corresponding to the size of the void. 
Therefore, by measuring the normalized shear, future galaxy surveys will be able not 
only to determine whether we live in a homogeneous universe or not, but to measure 
the size and depth of the void. 

In section 2 we describe the general LTB void models, giving the corresponding 
Einstein-Friedmann equations, as well as parameterizations of their solutions. In a 
subsection we describe the GBH constrained model, where we assume the Big Bang is 
homogenous, and thus the model only depends on a single function, the inhomogeneous 
matter ratio QM{f)- In section 3 we introduce a new observable, the normalized cosmic 
shear, which could help distinguish observationally homogeneous FRW models from 
arbitrary inhomogeneous universes. In section 4 we show how to properly calculate 
the size of standard rulers in LTB models, constructing an effective Einstein de 
Sitter background, that can be used for calculating early universe quantities in the 
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corresponding LTB model using standard formulae. In section 5 we analyse present 
observations of BAO along the line of sight, as well as the recent UNION compilation 
of Supernovae, and give constraints on the model from current observations. Finally, in 
section 6 we give a discussion of future prospects and some conclusions. 



2. Lemaitre-Tolman-Bondi void models 



The LTB model describes general radially symmetric space-times and can be used as a 
toy model for describing voids in the universe [38l [391 SO] • The metric is 

ds^ = -dt^ + X\r, t) dr^ + A\r, t) dQ^ , (1) 

where dQ"^ = dO"^ + sm^Odcj)^. Assuming a spherically symmetric matter source with 
negligible pressure, 

Tt = -pM(r,t)5o^5°, (2) 

the (0,r) component of the Einstein equations, G^^ = 0, implies X{r,t) = 
A'(r, t)/ ^Jl — k(r), with an arbitrary function k(r) playing the role of the spatial 
curvature parameter. The other components of the Einstein equations read flEi [T9l [23] 

+ 2HtHl + ^ + ^ = Svr G pm , (3) 

2Ht + ?>HI + ^ = Q, (4) 

where = dt ^ ' = dr, and we have defined the transverse and longitudinal Hubble rates 
as Ht = A/A, and Hl = A' /A'. Integrating the last equation, we get 
2 _ F(r) k{r) 

with another arbitrary function F{r), playing the role of effective matter content, which 
substituted into the first equation gives 
F'(r) 



871 GpM{r,t), (6) 



A'A^{r,t) 

where pM{r,t) is the physical matter density. 

The boundary condition functions F{r) and k{r) are specified by the nature of the 
inhomogeneities through the local Hubble rate, the local total energy density and the 
local spatial curvature, 

F{r) = H',{r)QM{r)Al{r), (7) 

kir) = H',ir)(nMir)-l)Al{r), (8) 



where functions with subscripts correspond to present day values, Ao(r) = yl(r, to) 
and Ho{r) = HTir^to). With these definitions, the r-dependent transversal Hubble rate 
can be written as [iHl [I9] 



H^{r,t) = H',{r 



^A/(r) . + (1 - UM{r)) 



A{r,t) J ' '''' " \A{r,t) 



(9) 
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and we fix the gauge by setting v4o(r) = r. For fixed r the above equation is equivalent 
to the Friedmann equation of an open universe, and have an exact parametric solution, 
while also very good approximate solutions can be found by Taylor expanding around 
an Einstein de Sitter solution (see [23] for details). 

For light travelling along radial null geodesies, ds"^ = dVL^ = 0, we have 

? = (10) 
dr ^1 - k{r) 

which, together with the redshift equation, 

cilog(l + ^) _ ^ A'{r,t) 
dr ^1 - k{r) 

can be written as a parametric set of differential equations, with = log(l + z) being 
the effective number of e-folds before the present time, 
dt A'{r,t) 



dN A'{r,t) 



(12) 



dr_ ^ ^ ^l-k{r) 
dN A'{r,t) 

Notice, that while the angular diameter distance dA = A[r,t) only depends on the 
integral of the transversal Hubble rate HT{r,t), the comoving and the luminosity 
distances di = {1 + z) = (1 + z^dA depend through Eq. (|T2|) on a mixture of 
Ht and Hi. 



2.1. The constrained CBH model 

In general LTB models are uniquely specified by the two functions k{r) and F{r) or 
equivalently by Ho{r) and f2jv/(r), but to test them against data we have to parameterise 
the functions, to reduce the degrees of freedom to a discrete set of parameters. For 
simplicity in this paper we will use the constrained GBH model [23] to describe the 
void profile. First of all, it uses a minimum set of parameters to make a realistic 
void profile, and secondly, by construction, the time to the Big Bang is constant for 
spatial slices, as would be expected in a generic model of inflation. Moreover, void 
models with an inhomogeneous Big Bang contain a mixture of growing and decaying 
modes, and consequently the void does not disappear at high redshift making them 
incompatible with the Standard Big Bang scenario [H]. If we only consider constrained 
LTB models, then at high redshifts and/or at large distances the central void is reduced 
to an insignificant perturbation in an otherwise homogeneous universe described by an 
FRW metric, and physical results for the early universe derived for FRW space-times 
still hold, even though we are considering an LTB space-time. 

The second condition gives a relation between Ho[r) and QM{f), and hence 
constrain the models to one free function, and a proportionality constant describing 
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the overall expansion rate. Our chosen model is thus given by 




(14) 







1 



sinh 



n=0 



(2n+ l)(2n + 3) ' 



2[fi/<(r)] 



(15) 



where flxif) = 1 — ^ni^), and the second equation follows from the requirement of a 
constant time to a homogeneous Big Bang. We use an "inflationary prior" , and assume 
that space is asymptotically fiat, i.e. in the following we set Qont = 1- The model has 
then only four free parameters: The overall expansion rate Hq, the underdensity at the 
centre of the void f2in, the size of the void ro, and the transition width of the void profile 
Ar. For more details on the model see Ref. [23] . 

3. The Raychaudhuri equation and the normalized shear 

Now that cosmological surveys are beginning to provide detailed maps of the universe up 
to significant distances we can start to consider measuring deviations from homogeneous 
FRW models. The main difference between FRW and LTB models is that the latter 
introduces two different components to the rate of expansion: There is both a transverse 
(perpendicular to the line of sight) and longitudinal (along the line of sight) components. 
This induces a differential growth of the local volume of the universe. Let us see how 
to quantify this geometrical distortion in the expansion of the universe. 

Suppose we start with a set of comoving observers localized on a sphere of volume 
in a LTB universe. Each observer has a unit vector ra^ tangent to its trajectory, 
and follows a geodesic of the metric, i.e. n^n".^ = 0. If we now consider a congruence 
of observers (i.e. geodesies) we can follow the change in positions of the observers 
with time. The covariant derivative of the congruence can be split into three distinct 
components. 



where P'^ = 6'^^+n^n,y projects out any tensor into a plane orthogonal to the congruence 
n'^. The tensor Q^^, measures the extent to which neighbouring trajectories deviate 
from remaining parallel. The three components of Q^i, in Eq. ( |T6l) have different 
physical meanings. The trace gives the expansion rate of the congruence, G = n'^.^, 
and characterizes the growth of the overall volume V of the sphere of observers. The 
traceless symmetric part is the shear tensor of the congruence. 



and represents the distorsion in the shape of the sphere of test particles (observers), 
giving rise to an ellipsoid with possibly different axes' lengths but same volume V. 
Finally, the antisymmetric part is the vorticity tensor. 




(16) 




(17) 
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which describes the rotation of the sphere of test particles around the center, again 
without changing the volume of the sphere. Note that these are purely geometrical 
quantities and are independent of the dynamics responsible for time evolution in this 
metric. By computing the covariant derivative of B^,y along the worldlines of observers 
and taking the trace, one finds the famous Raychaudhuri equation, 

^ = -^e^ - a.^a^" + u.^uj'^'' - R.^n^n" . (19) 

The last term of the equation is the only one sensitive to the theory of gravity. 
For general relativity and geodesic observers it becomes R^yU^n^ = 'SirGiT^y — 
}^g^yT)n^n'^ = 47rG(p + 3p). It is this term which allows for accelerated expansion 
in a ACDM-FRW model (for which both shear and vorticity vanish). 

We can now compute all the terms in this equation in our LTB models. The 
global rate of expansion becomes = Hi + 2Ht', the spatial shear is non-vanishing, 
aij = {Ht — Hi) Mij, where Mj = diag(— 2/3, 1/3, 1/3) is a traceless symmetric matrix; 
while the vorticity tensor vanishes, = 0. Note that the shear component of the 
congruence vanishes asymptotically as we approach the Einstein-de Sitter universe, as 
well as locally, close to the center of the LTB void. However at intermediate redshifts 
it has its maximal deviation from zero, where the longitudinal and transverse rates of 
expansion differ maximally, corresponding to the edge of the void. This shear component 
is a distinctive feature of inhomogeneous cosmological models, in fact it is one that could 
be used to distinguish FRW from LTB models in the near future with the next generation 
of astronomical surveys like DES, PAU, BOSS and LSST. 

In order to be quantitative, we should find a directly measurable quantity that can 
be computed and compared with observations. From a theoretical point of view, the 
natural variable to consider is the ratio of shear to expansion, 

where = Cija^^ = ^{Ht — . However, these rates are difficult to measure directly 
(in particular the transverse rate of expansion Ht can only be measured indirectly). 
More appropriate would be to use physical quantities like angular diameter distances 
and radial BAO scales, as a function of redshift or its derivatives. 

One can indeed find a relation between the variable e and those quantities, by 
noting that in LTB models, the angular diameter distance dA{z) is nothing but A{r,t) 
in the metric ([I]). Its derivatives give the various rates of expansion, Ht = A/A and 
Hl = A' /A', but note that in general, the time derivatives in an inhomogeneous universe 
are not directly related to the redshift derivatives since there is also a contribution from 
spatial gradients, 

d _ v^l - k{r) d 1 d 
In ~ A'Hl d^ ~ Ihdi' 
where we have used Eqs. f ll2tfT^ . and therefore H^ and Hx are not directly measurable 
quantities, contrary to the case in FRW models. Computing explicitly the derivatives 
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of A{r,t) with respect to redshift, we find 

Ht _ ^l^W) d\nA _ y/l - k{r) dlnd^(^) 

Hl AHl dN dAiz)HL{z) dln(l + z)' ^ ' 

which gives for the new variable, the normahzed shear e, as a function of the hne of 

sight rate of expansion and the angular diameter distance, 

^ ^ VT^^ -HL{z)d,[il + z)dA{z)] 

3HLiz)dAiz) + 2^1 - k{r) - 2Hl{z) + z) c/^l^)] ' 
This quantity still depends on an unknown, the spatial curvature k{r). However, we 
can integrate Eq. ( fT3l) to give 

1 dz' 

which is not far from unity in a wide range of LTB models that agree with observations, 
see Ref. [23| 12^. In this case, the normalized shear becomes 
, . ^ l-H,{z)d,[{l + z)dA{z)] 

^ ' ^ mL{z)dA{z) + 2 - 2Hl{z) + z) dA{z)] ' ^ ^ 

which can in principle be measured by any of the future surveys. In particular, the 
rate of expansion along the line of sight has recently been measured by Gaztanaga et 
al. at two different redshifts using the radial BAO scale EZ], and together with 
measurements of angular diameter distances could give a value for e{z). Note that in 
FRW universes the normalized shear ([22]) vanishes identically since = Hx = H and 




(1 + z)dA{z) = sinh 

-noV "A- 

Note also that the function HL{z)dA{z) in the denominator is nothing but the Alcock- 
Paczynski factor f{z), which is used as a geometric test for the existence of vacuum 
energy in ACDM FRW models. On the other hand, it is worth making emphasis on 
the fact that this normalized shear ( l23l) is independent of the value of Ho, still a mayor 
uncertainty in cosmology. 

In Fig. [1] we have plotted the normalized shear e for a LTB model with cosmo logical 
parameters corresponding to the best fit model that accommodates current CMB, SNIa 
and RBAO observations (see section [5]). A survey that can measure e in Eq. (|23ll 
significantly above zero will be able to distinguish homogeneous FRW models from 
inhomogeneous universes. The maximum value of the shear is attained when the density 
gradient is maximal, and thus the redshift of the maximum gives information about the 
size of the void. Note that the typical value of the normalized shear can be small, of 
order a few percent for LTB models, so this may impose severe requirements on future 
proposed surveys. A way to characterize the ability of a survey to distinguish between 
homogenous and inhomogeneous models of the Universe would be to ask them to detect 
a nonvanishing shear (e.g. max(£) ~ 0.03) at 95% confidence level, see inset of Fig. [H 
Note also that the sign of the (normalized) shear parameter can be used to characterize 
and distinguish a void from an overdensity. 
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Figure 1. The redshift dependenee of the normalized slrear ([55]) for an LTB universe 
with cosmological parameters eorresponding to the best-fit model (full line) and to the 
model with minimum (dashed line) with 1. 2, and 3-cr envelopes, that accommodates 
current CMB, SNIa and RBAO observations. We also plot in the inset the distribution 
function of maximum shear in 1, 1+2 and 1+2+3-ct bins. 



4. Standard Rulers in LTB models 

Many bounds from observational cosmology, such as the sound horizon (a "standard 
ruler"), the CMB, and applications of the Alcock-Paczinsky test, are derived by 
considering and calculating physical scales and processes in the early universe, and 
are based on the implicit assumption of an underlying FRW metric. To test LTB 
models against these observational data we have to connect distance scales, redshifts, 
and expansion rates in the early universe to those observed today. 

Starting from an approximately uniform universe at a high redshift z^. in the LTB 
model, the expansion rate and matter density becomes gradually inhomogeneous, and 
a uniform co moving physical length scale / in the early universe at Ze-, for example the 
sound horizon, is not uniform at some later redshift z, because of the space dependent 
expansion rate. In particular the comoving length at the current time to depends on 
how much relative expansion there has been at different positions since the formation 
of the uniform length scale 



where r^ is the radial coordinate of an observer far away from the void. While it is clear 
that the physical length becomes scale dependent due to the inhomogeneous expansion 
of the Universe, the scale-dependence of the comoving scale is a consequence of it being 
measured at to- If instead we defined the comoving length scales to be measured in 
the early universe at t{ze), then indeed l{r{z)) would be independent of the observer 




(24) 
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position. The convenience of the above formula is, that the LTB models we consider in 
this paper are asymptotically Einstein de Sitter, and we can easily compute comoving 
scales at infinity. 

Not only do we have to take into account that scales are different, but we would 
also like to use the standard framework to compute early universe quantities. Hence, 
we need to compute the equivalent Hubble constant at the current redshift, that an 
observer in a Einstein de Sitter universe would have now at to to have the same redshift 
history far away from the void, that the central observer has in the LTB model. Using 
the normal FRW equation for an Einstein de Sitter universe we can write: 

i/2(^e) = ^eff(l + ^ef (25) 

where Ze is a high enough redshift where the lightcone is far away from the void, but 
where the radiation density has not become significant yet, typically z ~ 100. The 
l.h.s. has to be evaluated in the LTB model, while the r.h.s. gives the equivalent Einstein 
de Sitter Hubble constant H^s, that can be used as input for calculating e.g. the sound 
horizon using standard methodi§|. 



5. Constraints from observations 



To constrain the parameters of our model we use three data sets: The first acoustic peak 
in the CMB as measured by WMAP |12], the Union Supernovae data set [S] consisting 
of 307 Type la supernovae, and the radial baryon acoustic scale measured at z = 0.24 
and z = 0.43 using large red galaxies in the SDSS DR6 [SUl EZ]- 

The likelihood for the acoustic CMB scale is calculated as in [23]. For the 
UNION supernovae we have used the full covariance matrix published by the Supernova 
Cosmology Project [HI US] including systematic errors, and we have adapted their 
likelihood code to our easyLTB program. We find the radial baryon acoustic differential 
redshift Az in our model as 

A.,,, = MlllM (26) 

c 

and compare it to the measurements, including systematic errors, published in [37] . The 
size of the sound horizon is redshift dependent and computed as detailed in section HI 

Several changes can be observed (see Fig. HI) compared to our previous analysis in 
[23] : Clearly the UNION supernovae prefer a higher overall density contrast of the void, 
compared to the Davis et al. compilation [H]. This is curious, given that most of the 
supernovae are overlapping between the two data sets, and it is not a priori clear where 
the difference is. A possible cause, apart from the increased number of supernovae in the 
UNION set, could be that all supernovae in the UNION compilation has been reanalysed 
with the SALT light curve fitter, presumably giving a more homogeneous data set 

§ See |23j for the general case when the space time is not asymptotically Einstein de Sitter, and Ref. |25j 
for a treatment, that also fixes the angular diameter distance to the correct value in the effective Einstein 
de Sitter model, so the full angular powerspectrum can be computed by standard CMB packages 
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Hq Hr=0 Hr=oo 


^in 


ro Ar 


units 


100 km Mpc^i 




Gpc ro 


Priors 


0.50-0.95 0.4-0.89 0.33-0.63 


0.05-0.35 


0.5-4.5 0.1-0.9 


Best Fit ± 2-a 


0.67±0.03 0.58 0.45 


0.16±0.09 


2.7±0.8 0.44(>0.12) 


Minimum 


0.68 0.59 0.45 


0.145 


2.35 0.85 



Table 1. Priors used when scanning the parameters of the model, the marginaUsed 
best fit values with 2-a confidence levels (see Fig. [5]), and the values at the minimum 
X^. Hq is only a pre factor for HQ{r) and the priors on Hr=o and iJr=oo are derived 
from the priors on fli^ and Hq- Because they have a complicated prior, we give them 
without confidence limits. 
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Redshift 

Figure 2. The 1,2, and 3-<t limits for the sound horizon differential redshift Az along 
the line of sight, the longitudinal Hubble rate Hl, and the sound horizon together with 
the best fit model (full line). The observed Az with errors are given as boxes, and the 
equivalent best fit ACDM model from is indicated by a dashed line. 

without spurious offsets between different groups of supernovae. Tlie magnitude-redsliift 
relation is sliown in Fig. [3l and even tliougli we only imposed as good a fit as possible 
to the data, the curves of the LTB models with the best fit, and the minimum are 
both almost identical to the best fit ACDM model, illuminating the very small freedom 
in the magnitude-redshift relation allowed by current data. Notice also that at z > 2 
the slopes of the three curves are approximately the same, because at those redshifts all 
models are effectively Einstein de Sitter. Hence, a ~5 Gpc sized void would be able to 
mimic the ACDM magnitude-redshift relation up to any redshift. 

The radial baryon acoustic (RBAO) scale is pushing the void size upwards compared 
to the isotropised BAO data of Percival et al [13]. This can be understood in terms of 
the higher redshift of the second data point (z=0.45), compared to the second data point 
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Figure 3. The UNION supernova data together with the magnitude-redshift relations 
for the best fit and minimum LTB models compared to the standard ACDM model. 

of Percival et al (z=0.35). Given that the RBAO measurements are in perfect agreement 
with a ACDM model, the void has to be bigger, to have the same RBAO differential 
redshifts, Az, at the two different redshifts. In Fig. [2] is shown the envelope curves of 
A2;, rs{z), and Hl{z) as a function of redshift derived from models which are inside the 
1, 2, and 3-0" likelihood contours. While the best-fit LTB model is perfectly inside the 
error bars of the two data points, it is very different from the ACDM model at higher 
and lower redshifts, and just another redshift bin, outside current measurements, would 
help considerably in discriminating between the models. While Az is monotonically 
increasing in ACDM, it is almost constant a.t z < 0.5 in the best-fit LTB model. This 
can be understood by looking at the inserts in Fig. O The sound horizon is much larger 
at lower redshift, because of the higher expansion rate near the centre of the void, 
and this is almost compensated by the longitudinal Hubble rate, giving a practically 
fiat RBAO differential redshift. Unfortunately it is difficult due to measure the RBAO 
feature at lower redshifts, because of cosmic variance, but in the near future dedicated 
RBAO surveys will detect it at higher redshifts. 

The CMB acoustic scale fixes the overall Hubble parameter Hq, because our 
models are asymptotically flat, and because we have fixed cubaryon = 0.0223 to the 
best fit WMAP3 value. Overall the model is capable of yielding a very good fit 
to data with a minimum = 307.7 for 310-4 d.o.f. In figure O we show the one 
dimensional likelihoods, and it can be seen that the different data sets yield consistent 
and overlapping marginalised likelihoods, further supporting that the model is a very 
good and consistent fit to the observations. 

6. Discussion and conclusions 

As an alternative explanation of the dimming of distant supernovae it has recently been 
advocated that we live in a special place in the Universe near the centre of a large void. 
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Figure 4. Likelihood contours for the combined data set is shown in yeUow with 1-, 
2-, and 3-cr contours, while the individual likelihoods for the SNIa, BAO, and CMB 
data sets are shown in blue, purple, and green respectively with 1- and 2-ct contours. 
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Figure 5. Marginalised one dimensional likcliliood curves for the combined data set 
is shown in yellow (full line), while the individual likelihoods for the SNIa, BAO, and 
CMB data sets are shown in blue (dot-dashed), purple (dot-dot-dashed), and green 
(dashed) respectively. The dotted lines indicate the 1-, 2-, and 3-(t levels, while the 
long dashed lines gives the best fit and the 2-a values. 

In this scenario, the universe is no longer homogeneous and isotropic and the apparent 
late time acceleration is actually a consequence of spatial gradients in the metric. 

If, at the end, local spatial curvature explains away the need for a cosmological 
constant we would be living through a paradigm shift similar to the one that occurred 
with the discovery of the expansion of the universe by Hubble and others. At that 
time, the metric of the universe was assumed to have maximal symmetry (time and 
space translations plus rotations), suggesting a static, homogeneous and isotropic 
universe, in agreement with the Perfect Cosmological Principle. The observed expansion 
of the universe removed the need for time translation invariance and only spatial 
sections were maximally symmetric, and we are left with a homogeneous and isotropic 
universe satisfying the usual Cosmological Principle, which states that any observer 
in the universe is equivalent to any other one. However, if the universe around us is 
actually curved, within a large void in an otherwise EdS universe, then it is no longer 
homogeneous, and we are left only with rotational invariance. The metric corresponds to 
that of a Lemaitre-Tolman-Bondi model. These models can accommodate the observed 
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dimming of distant supernovae, without the need for a cosmological constant, because 
spatial curvature also increases the luminosity distance. In fact, LTB models fit well 
the present SN data, even better than ACDM for some models [26] . 

In this paper we have tested a class of simple LTB models, against different 
cosmological data, that probes the overall geometry of the universe, as observed in 
the lightcone, namely: The Union supernovae data set of 307 Type la SNe [5], the 
recent measurement of the radial B AO scale in two different redshift bins [SU] , and the 
CMB acoustic scale as measured by WMAP [12]. We have shown that LTB models yield 
a convincing fit to current observations with only 4 parameters in the model, and an 
excellent for ths best-fit model of x^/d.o.f. = 307.7/(310 — 4). While not ruling out 
the void model, the observations do require a void of at least > 2 Gpc (at 95% C.L.), 
and a central underdensity that is in accordance with local galaxy survey measurements 
of fiin ~ 0.16 — 0.25. The predicted local Hubble parameter of 59 km s~^ Mpc~^ is 
within 2-cr of the Hubble key project value of 72 ± 8 km s~^ Mpc~^. 

Moreover, there are new specific predictions that these inhomogeneous models give 
for physical observables, the most important one being that the overall rate of expansion 
is not homogeneous and thus longitudinal and transverse expansion rates along the 
line of sight are different. In the language of differential geometry, this implies that 
the congruence of geodesic observers has a non-vanishing shear. If we normalize this 
shear w.r.t. to the expansion (they both appear on equal footing in the Raychaudhuri 
equation that governs the acceleration of the universe) then the normalized shear can 
be used as an observable for future surveys. We have found a very compact form of 
the normalized shear that can be measured in principle, written in terms of directly 
observable quantities. In fact, if future surveys measure at high confidence level a 
non-vanishing value for this new observable, then we will have to accept that we live 
in an inhomogeneous universe, since any homogeneous model (be it ACDM or any 
other generalized vacuum energy models) has vanishing shear. In a particular class of 
inhomogeneous models describing large spherical voids, also known as LTB models, the 
normalized shear has a maximum at the position of the steepest gradient in the void 
profile. Supernovae observations are compatible with an LTB model with void size of 
order a few gigaparsecs, well beyond present galaxy catalogs but within reach of the 
next generation. In those models, the maximum shear is of order a few percent, and 
may be difficult to measure over the whole redshift range. In the future, the ability of a 
galaxy survey to measure the normalised shear could be used as a discriminator between 
different proposals, in a way similar to (and independent of) the figure of merit for the 
characterization of the nature of dark energy in terms of an equation of state parameter 
and its derivative, as advocated by the Dark Energy Task Force. 

Remote measurements of the CMB anisotropies through precise determinations of 
the kinematic Sunyaev-Zel'dovich effect in a handful of clusters, as recently proposed by 
us [21], or via spectral distortions in the local CMB radiation induced by CMB photons 
passing through the void that are later rescattered back towards us during reionisation, 
proposed by Caldwell and Stebbins [16] are other new ways of testing the Copernican 



The radial BAO scale and Cosmic Shear, a new observable for Inhomogeneous CosmologieslQ 



Principle and constraining general void models. 

It is worth mentioning here the recent independent observations of a systematically 
large bulk flow on large scales both above and below 100 Mpc [171 HHl HH], based on 
the kinematic Sunyaev-Zel'dovich effect in X-ray clusters, and various peculiar velocity 
surveys in the IRAS-PSCz density field respectively. The size of the flow is significantly 
above that predicted in the Standard Model and constitutes a challenge for ACDM. 
Whether this new observation is due to a large void at distances of order gigaparsecs is 
still very uncertain, and we will probably have to wait for the next generation of galaxy 
surveys to rule out this possibility. 

The combination of the steadily increasing body of Type la Supernovae, 
constraining the transverse Hubble rate, together with upcoming galaxy surveys that 
will measure the radial BAO scale in many redshift slices, such as PAU and BOSS, 
constraining the longitudinal Hubble rate, will put significant pressure on LTB models, 
or measure for the first time the shear component of cosmic expansion. If not confirmed 
by observations, then only a small corner of parameter space - with very smooth void 
profiles and Ht ^ - will be left. 

Note added: When we were finalizing this work, we learned of a similar analysis 
made by Zibin et al. [25], where they also used the UNION supernova and the radial 
BAO of Ref. [Sni EZ] to constrain void models. They also tested their void model against 
the CMB anisotropies, and used a toy model for the ISW effect. They get similar void 
sizes, albeit they have bigger difficulties in fitting the data, due to their inclusion of the 
anisotropies in the CMB. 
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